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We consider a minimal formulation of SO(10) Grand Unified Theory wherein all the 
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mechanism, the large v^ — Vr mixing as a consequence of b— r unification at the GUT scale. In 

in ■ 

this picture, however, the CKM phase S lies preferentially in the second quadrant, in contra 



fermion masses arise from Yukawa couplings involving one 126 and one 10 of Higgs multi- 
plets. It has recently been recognized that such theories can explain, via the type-II seesaw 



diction with experimental measurements. We revisit this minimal model and show that the 
conventional type-I seesaw mechanism generates phenomenologically viable neutrino masses 

Ph. 

and mixings, while being consistent with CKM CP violation. We also present improved fits 



in the type-II seesaw scenario and suggest fully consistent fits in a mixed scenario. 



I. INTRODUCTION 

Grand Unified Theories (GUT) provide a natural framework to understand the properties of 
fundamental particles such as their charges and masses. GUT models based on SO(W) gauge 
symmetry have a number of particularly appealing features. All the fermions in a family fit in a 
single 16-dimensional spinor multiplet of SO(10). In order to complete this multiplet, a right- 
handed neutrino field is required, which would pave the way for the seesaw mechanism which 
explains the smallness of left-handed neutrino masses. 50(10) contains SU(5) and the left-right 
symmetric Pati-Salam symmetry group SU(4)c x SU{2)l x SU{2)h as subgroups, both with 
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very interesting properties from a phenomenological perspective. With low energy supersymmetry, 
50(10) and SU(5) models also lead remarkably to the unification of the three Standard Model 
gauge couplings at a scale Mqut ~ 2 x 10 16 GeV. 

In grand unified theories, the gauge sector and the fermionic matter sector are generally quite 
simple. However, the same is not true of the Higgs sector. Since the larger symmetry needs to 
be broken down to the Standard Model, generally one needs to introduce a large number of Higgs 
multiplets, with different symmetry properties under gauge transformations. If all of these Higgs 
fields couple to the fermion sector, one would lose much of the predictive power of the theory in 
the masses and mixings of quarks and leptons, and so also one of the attractive aspects of GUTs. 

Of interest then are the so-called minimal unification theories, in which only a small number 
of Higgs multiplets couple to the fermionic sector. One such realization is the minimal 50(10) 
GUT []| in which only one 10 and one 126 of Higgs fields couple to the fermions. These two Higgs 
fields are responsible for giving masses to all the fermions of the theory, including large Majorana 
masses to the right-handed neutrinos. This model is minimal in the following sense. The fermions 
belong to the 16 of 50(10), and the fermion bilinears are given by 16 x 16 = 10 s + 120 a + 126 s . 
Thus 10, 120 and 126 Higgs fields can have renormalizable Yukawa couplings. If only one of these 
Higgs fields is employed, there would be no family mixings, so two is the minimal set. The 126 
has certain advantages. It contains a Standard Model singlet field and so can break 5O(10) down 
to SU(5), changing the rank of the group. Its Yukawa couplings to the fermions also provide large 
Majorana masses to the right-handed neutrinos leading to the seesaw mechanism. It was noted 
in Ref. Q] that due to the cross couplings between the 126 and the 10 Higgs fields, the Standard 
Model doublet fields contained in the 126 will acquire vacuum expectation values (VEVs) along 
with the VEVs of the Higgs doublets from the 10. The 126 Yukawa coupling matrix will then 
contribute both to the Dirac masses of quarks and leptons, as well as to the Majorana masses of 
the right-handed neutrinos. 

It is not difficult to realize that this minimal model is highly constrained in explaining the 
fermion masses and mixings. There are two complex symmetric Yukawa coupling matrices, one of 
which can be taken to be real and diagonal without loss of generality. These matrices have 9 real 
parameters and six phases. The mass matrices also depend on two ratios of VEVs, leading to 11 
magnitudes and six phases in total in the quark and lepton sector, to be compared with the 13 
observables (9 masses, 3 CKM mixings and one CP phase). Since the phases are constrained to be 
between —1 and +1, this system does provide restrictions. More importantly, once a fit is found 
for the charged fermions, the neutrino sector is fixed in this model. It is not obvious at all that 
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including the neutrino sector the model can be phenomenologically viable. 

Early analyses found that just fitting the lepton-quark sector is highly constraining. Also, 

this fitting has been found to be highly nontrivial (in terms of complexity); therefore these analyses 
were done in the limit when the phases involved are either zero or ir. In such a framework, one finds 
that the parameters of the models are more or less determined by the fit to the lepton-quark sector 
(the quark masses themselves are not known with great precision, so there is still some room for 
small variations of the parameters). As a consequence, one could more or less predict the neutrino 
masses and mixings; however, since neutrino data was rather scarce at the time, one could not 
impose meaningful constraints on the minimal 50 ( 10) model from these predictions. 

In view of the new information on the neutrino sector gathered in the past few years Q, U, Q] , 
one should ask if this model is still consistent with experimental data. Interest in the study of this 
model has also been reawakened by the observation that b — r unification at GUT scale implies 
large (even close to maximal) mixing in the 2-3 sector of the neutrino mass matrix 6], provided 
that the dominant contribution to the neutrino mass is from type-II seesaw. There has been a 
number of recent papers studying the minimal 50(10) using varying approaches: some analytical, 
concentrating on the 23 neutrino sector [(], Q] , some numerical, either in the approximation that 
the phases involved in reconstructing the lepton sector are zero 
into account 



, or taking these phases 



L2j, 3- ons of these analyses seem to be that the minimal 50(10) cannot 

account by itself for the observed neutrino sector (although it comes pretty close). However, one 
might restore agreement with the neutrino data if one slightly modifies the minimal 50(10); for 
example, one can set the quark sector CKM phase to lie in the second quadrant, and rely on 
new contributions from the SUSY breaking sector in order to explain data on quark CP violation 
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131 ] . or even another Higgs 
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[l^ : or one might add higher dimensional opovalors lo the theory 
multiplet (a 120) which will serve as small perturbation to fermion masses 

In this paper we propose to revisit the analysis for the minimal 5O(10) model, with no extra 
fields added. The argument for this endeavor is that our approach is different in two significant 
ways from previous analyses. First, we use a different method than 3 13 i n fitting for the 
lepton-quark sector. Since this fit is technically rather difficult, and moreover, since the results 
of this fit define the parameter space in which one can search for an acceptable prediction for the 
neutrino sector, we think that it is important to have an alternative approach. Second, rather 
than relying on precomputed values of quark sector parameters at GUT scale, we use as inputs 
Mz scale values, and run them up to unification scale. This allows for more flexibility and we 
think more reliable predictions for the parameter values at GUT scale. With these modifications 







in our approach, we find that we agree with some results obtained in [12 . 1 1 31 ] (in particular, the 
fact that type-II seesaw does not work well when the CKM phase is in the first quadrant), but not 
with others. Most interesting, we find that it is possible to fit the neutrino sector in the minimal 
50(10) model, in the case when type-I seesaw contribution to neutrino mass dominates. We also 
present a mixed scenario which gives excellent agreement with the neutrino data. 

The paper is organized as follows. In the next section we give a quick overview of the features 
of the minimal 50(10) model relevant for our purpose. In section III we address the problem of 
fitting the lepton-quark sector in this framework. We also define the experimentally allowed range 
in which the input parameters (quark and lepton masses at Mz scale) are allowed to vary. We 
start section IV with a quick overview of the phenomenological constraints on the neutrino sector. 
There we provide a very good fit to all the fermion masses and mixings using type-I seesaw. We 
follow by analyzing the predictions of the minimal 50(10) model in the case when type-II seesaw is 
the dominant contribution to neutrino masses. We then analyze the predictions in a type-I seesaw 
dominance scenario, and in a scenario when both contributions (type-I and type-II) have roughly 
the same magnitude. We end with our conclusions in Sec. V. 



II. THE MINIMAL SO(10) MODEL 

The model we consider in this paper is an super symmetric 50(10) model where the masses of the 
fermions are given by coupling with only two Higgs multiplets: a 10 and a 126 Q]. Both the Hio 
and Hj26 contain Higgs multiplets which are (2,2) under the SU(2)l x SU(2)r subgroup. Most 
of these (2,2) Higgses acquire mass at the GUT scale. However, one pair of Higgs doublets H u and 
Hd (which generally are linear combinations of the original ones) will stay light. (Details about the 
Higgs multiplet decomposition and 50(10) breaking can be found, for example, in 17, 18, 19, 2o|). 
Upon breaking of the SU(2)l X U(l)y symmetry of the Standard Model, the vacuum expectation 
value of the H u doublet will give mass to the up-type quarks and will generate a Dirac mass term 
for the neutrinos, while the vacuum expectation value of the doublet will give mass to the 
down-type quarks and the charged leptons. 

The mass matrices for quarks and leptons wil then have the following form: 

M u = KyYlQ + «4Yi 2 6 

M d = KdYm + K d Y 126 
M? = k u Y 10 - 3k' u Y 126 
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Ml = K d Y W - 3/^*126 (1) 

where lio,^i26 are the Yukawa coefficients for the coupling of the fermions to the Hio and Hj^g 
multiplets respectively. Note that in the above equations the parameters n Ut d, k' u d as well as the 
Yukawa matrices are in general complex, thus insuring that the fermion mass matrices will contain 
CP violating phases. 

The SO(10) Hj^g multiplet also contains a (10,3,1) and a (10,1,3) Pati-Salam multiplets. 
The Higgs fields which are color singlets and SU(2)r/ SU{2)i triplets (denoted by A# and A^) 
may provide Majorana mass term for the right-handed and the left-handed neutrinos. One then 
has: 

M uR = (A R )Y 126 

M vL = (A L ) Y 126 . (2) 

If the vacuum expectation of the Ar triplet is around 10 14 GeV then the Majorana mass term for 
the right-handed neutrinos will give rise, through the seesaw mechanism, to left-handed neutrino 
masses of order eV. On the other hand, the VEV of A^ contributes directly to the left-handed 
neutrino mass matrix (this contribution is called type-II seesaw), so this requires that the (A/,) is 
either zero or at most of order eV. This requirement is satisfied naturally in such models, since A^, 



generally aquires a VEV of order M|/(A^j) 21]. 



III. LEPTON AND QUARK MASSES AND MIXINGS 

Our first task is to account for the observed lepton and quark masses, and for the measured 
values of the CKM matrix elements. By expressing the Yukawa matrices Yio and Y126 in Eqs. Q 
in favor of M u and M^, we get a linear relation between the lepton and quark mass matrices; at 
GUT scale: 

Mi = aM u + bM d , (3) 

where a and b are a combinations of the K u ^x! ud parameters in Eq. For simplicity let's 

work in a basis where is diagonal (this can be done without loss of generality). Then, 
M n = Vc KM M^ Vckm with the diagonal up-quark mass matrix. If we allow the en- 
tries in the diagonal quark mass matrices to be complex: = diag{m u e tau ,m c e tac ,mte iat }, 
M d = diag{m,ie Lbd ,m s e lbs ,mi,e lbb }, then the CKM matrix can be written in its standard form as a 
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function of three real angles and a phase: 



Vckm 



1 C12C13 S12C13 si 3 e~ tS ^ 

-S12C23 - CnS23Sl3e lS C12C23 - Sl2S23Sl3e lS S23C13 
y S12S23 - Cl2C23Sl3e 4<5 -C12S23 ~ •Sl2C23'Sl3e*' 5 C23C13 J 



(4) 



Since their phases can be absorbed in the definitions of the Oj,6j parameters, we will take the 
coefficients a, b to be real, too. One of the quark mass phases can be set to zero without loss of 
generality, we set bb = 0. It should be noted that a common phase of {a, b}, which we denote as 
a, will appear in the Dirac and Majorana mass matrices of the neutrinos, and will be relevant to 
the study of neutrino oscillations. 

The relation (J2J) will generally impose some constraints on the masses of the quarks and leptons. 
For example, if we take all the phases to be zero (or tt), then on the right-hand side of the equation 
there are just two unknowns, the coefficients a, b. On the other hand, the eigenvalues of the lepton 
mass matrix are known, which will give us 3 equations. It is not obvious, then, that this system can 
be solved; however, early analysis [3, Q] shows that solutions exist, in the range of experimentally 
allowed values, provided that the quark masses satisfy some constraints. 

Newer studies |inl . 11, 3] allow for (some) phases to be non-zero, and thus relaxes somewhat 
the constraints on quark masses. However, it is interesting to note that these solutions are not 
very different from the purely real case. That is, most of the phases involved have to be close to 
zero (or tt), and the values of the parameters do not change by much. We shall explain this in the 
following. 

The algebraic problem of solving for the lepton masses in the case when the elements of the 
matrices M u , are complex is quite difficult. This would involve solving a system of 3 polynomial 
equations of degree six in unknown quantities a, b. Most of the analysis so far is done by numerical 
simulations (some analytical results are obtained for the case of 2nd and 3rd families only [fj,Q]). In 
this section we attempt to solve the full problem (with all the phases nonzero) in a semi-analytical 
manner, that is, by identifying the dominant terms in the equations and obtaining an approximate 
solution in the first step, which can then be made more accurate by successive iterations. 

Due to the hierarchy between the eigenvalues of the lepton mass matrix m e <C m„ <C m T one 
can suspect that the mass matrix itself has a hierarchical form. This assumption is supported by 
the observation that the off-diagonal elements of M; are indeed hierarchical; for example L13/L23 — 
V31/V31 <C 1. [Lij is a short-hand notation for (Mi)ij, Vij are the ij elements of Vckm-) Then, 
the three equations for the invariants of the real matrix LL^ (the trace, the determinant and the 



7 



sum of its 2 x 2 determinants) become: 

|L 3 3| 2 + 2|L 23 | 2 =i m 2 T 
\L 22 L 33 - L 2 23 \ 2 ~ mjm£ 

Det [LL f ] = m^m^mj . (5) 

We find it is convenient to work in terms of the dimensionless parameters rhi = rrii/m T (i = 
e,fj,),a = arh t ,b = brhb, and the ratios r c = m c /m t ,r s = m s /mb,r u = m u /m t ,rd = md/rub- 
Explicitly from the equations above in terms of these parameters, we obtain: 

L 33 = e iQl = aVl 3 + b e~ lZi 
K 23 = m^ = {br s + ar c e lz ^L 33 + dbVl 2 e^ bb - b ^ 
A = m e m^ = b r d A 23 - a 2 b e^~ b ^ (r s {V 3l V 33 ) 2 

+2r c V 31 V 32 V 21 V 22 e l{z ^ + r^V^V^e 1 ^ (6) 

Here we have kept only the leading terms, using Tui^d ^ciTs 1* Moreover, note that only 
phase differences like a, t — 6i, Oj — chj can be determinated from Eq. ©; therefore, by multiplying 
with overall phases, we have written Eqs. © in terms of these differences (with the notation 
Zi = ai - bi). 

The key to solving this system is to recognize that there is some tuning involved. Analyzing 
the first two equations leads to the conclusion that a,b — O(10). Then the phase z 3 in the first 
equation should be close to tt so that the two terms almost cancel each other. Similar cancellations 
happen in the second and the third equations, which require respectively that bb — b s ~ tt and 
bd — Of. Also, in the third equation, neglecting the small electron mass on the left hand side 
results in: 



a ^r s \V 3l V 33 \ 2 ■ [<) 



For values of the parameters in the experimentally feasible region, this is consistent with the above 
estimate a ~ 0(10). 

Analytically solving Eqs. (jHJ) with the approximations discussed will provide solutions for the 
phases and parameters a, b accurate to the 10% level. Using these first order results, one can 
compute and put the neglected terms back in Eqs. (|5ffi|) . which can be solved again, thus defining 



1 Note that taking these phase differences to tt or zero results in exactly the mass signs which the analysis in [f(l | 
found to work for the real masses case. 
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an iterative procedure which can be implemented numerically, and brings us arbitrarily close to 
the exact solution. We find that 5 to 10 iterations are usually sufficient to recover the fi and r 
masses with better than 0.1% accuracy (m e can be brought to a fixed value by multiplying with 
an overall coefficient). 

We end this section with some comments on the range of input parameters (masses and phases) 
which allow for a solution to Eq. @. As we discussed above, the phases are either close to ir or to 
zero. This is required by the necessity to almost cancel two large terms in the right-hand side of 
Eqs. ©. One can see that the larger the absolute magnitude of these terms (for example |aV| 3 | 
and b in the first equation), the more stringent are the constraints on the phases. The opposite 
is also true; the smaller the a and b parameters, the more the phases can deviate from ir, and 
generally the easier it is to solve the system. This means that lower values of 5, b are preferred; 
from Eq. 0, this implies a preference for low values of the ratio md/m s 2 (there is not much scope 
to vary m^). It turns out that low values of m s and large values of m c can also help, since they 
lower the absolute magnitude of the larger term on the right-hand side of the equation for A23 in 
©. Previous analysis found indeed that fitting for the lepton masses require a low value for m s 



A. Low scale values and RGE running 

As was discussed in the above section, the relation © implies some constraints on the quark 
masses (the lepton masses being taken as input). That is, not all values of quark masses consistent 
with the experimental results are also consistent with the model we use. Our purpose first is to 
identify these points in the parameter space defined by the experimentally alowed values for quark 
masses, 

Let us then define what this parameter space is. Altought the relations in the previous section 
hold at GUT scale, one must necessarily start with the low energy values for our parameters. 
We choose to use as input the values of the quark masses and the CKM angles at the Mz scale. 
Estimates of these quantities can be found for example in 22]. However, we consider some of 



their numbers rather too precise (for example, their error in estimating the masses of the s and 
c quarks are only 25%, respectively 15%, while the corresponding errors in PDG are much 
larger). Therefore, in the interest of making the parameter space as large as possible, we use the 



2 This also means higher values for | V31 1 are preferred. Since V31 = S12S23 — ci2S23Si3e 115 , this implies a preference 
for values of the CKM phase 8 close to ty (as noted in fT^ - ). 
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following values: 

• for the second family: 70 MeV < m s (M z ) < 95 MeV 3 ; 650 MeV < m c (M z ) < 850 MeV. 
With a running factor from Mz to 2 GeV of around 1.7, these limits would translate to 
values at 2 GeV scale of: 120 MeV < m s (2 GeV) < 160 MeV ; 1.1 GeV < m c (2 GeV) < 
1.44 GeV. Lattice estimations (m c /m s )2GeV — 12 [2J] would indicate a value in the lower 
part of the range for m s , and a upper part for m c . 

• for the light quarks: here generally the ratio of quark masses are more trustworthy than 
limits on the masses themselves; we therefore use use 17 < m s /m,i < 23 (as noted in the 
previous section, high values of this ratio are preferred), and 0.3 < m u /rrid < 0.7. We note 
here that m u is a parameter which does not affect the results much. 

• for the heavy quarks: 2.9 GeV < m b {M z ) < 3.11 GeV, (or 4.23 GeV < m b (m b ) < 4.54 
GeV) and for the pole top mass 171 GeV < Mi < 181 GeV (the corresponding MS mass is 
evaluated using the three loops relation, and comes out about 10 GeV smaller). 

• the CKM angles at Mz scale: 

sia = 0.222 ± 0.003 , s 23 = 0.04 ± 0.004 , s 13 = 0.0035 ± 0.0015. 

For the gauge coupling constants we take the following values at M z scale: a\{M z ) = 
1/58.97, a%{Mz) = 1/29.61, as (Mz) = 0.118. With these values at low scale one can get unifi- 
cation of coupling constants at the scale Mqut ~ 10 16 . The exact value of Mqut, as well as the 
values of the fermions Yukawas at the unification scale, will depend also on the supersymmetry 
breaking scale (Msusy) and tan/3, the ratio between the up-type and down-type SUSY Higgs 
VEVs. We generally consider values of Msusy between 200 GeV and 1 TeV, and tan/3 between 5 
and 60. 

Having chosen specific values of the parameters described above, we then run the fermion 
Yukawa coupling and the quark sector mixing angles, first from M z to Msusy, using two- loop 
Standard Model renormalization group equations; then we run from SUSY scale to the GUT scale 
using two loop 4 SUSY RGEs |26j. After computing the neutrino mass matrix at GUT scale, we 

3 Note that the lower limit for m s {Mz) is rather low compared with |2^|: however, the value at 2 GeV scale is well 
within the limits cited in [SJ. Lattice results also seem to favor smaller values of m s (2 GeV) [2e|. 

4 More precisely, we use the two-loop RGEs for the running of the gauge coupling constants and the third fam- 
ily fermions (6, t and r). To evaluate the light fermion masses, we use the one- loop equations for the ratios 
m Si i/mi, m CjU /m t and m M , e /m T . This approximation is justified, since the leading two-loop effect on the fermion 
masses comes from the change in the values of gauge coupling constants at two-loop; however, the contributions 
due to the gauge terms are family-independent and will not affect these ratios. 
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run its elements back to Z scale j^. and evaluate the resulting masses and mixing angles. 



IV. NEUTRINO MASSES AND MIXINGS 

In the present framework, there are two contributions to neutrino masses. First one has the 
canonical seesaw term: 

{M v ) seesawI = MfM^M? (8) 

with Mr = VrYx26 an d given by However, the existence in this model of the (10,3,1) 
Higgs multiplet implies the possibility of a direct left-handed neutrino mass term when the SU(2)l 
Higgs triplet A.l from this acquires a VEV vl (as it generally can be expected to happen). The 
neutrino mass contribution of such a term would be 

(M u ) seesawII = v L Y 126 = \M R (9) 



where vl = l v wea k/ v R anc ^ 7 is a factor depending on the specific form of the Higgs potential 21]. 

The scale of the canonical seesaw contribution Eq. (JSJ) (which we call type-I seesaw in the 
following) to the left handed neutrino mass matrix is given by v%j eak /vR. The contribution of the 
type- II seesaw factor (Eq. @) is of order lv^ eak jvR. One cannot know apriori how the factor 7 
compares with unity, therefore one cannot say which type of seesaw dominates (or if they are of 
the same order of magnitude). Therefore, in the following each case will be analyzed separately. 

However, let us first review the current experimental data on the neutrino mixing angles and 
mass splittings. Latest analysis j3] sets the following 3a bounds: 

• from ZAj — v T oscillations: 

1.4 x 10 -3 eV 2 < Am^ 3 < 3.3 x 10 -3 eV 2 ; 0.34 < sin 2 6» 23 < 0.66 ; 

with the best fit for Am^ = 2.2 x 10~ 3 eV 2 and sin 2 #23 = 0.5 (from atmospheric and K2K 
data). 

• from v e — Va oscillations: 

7.3 x 10" 5 eV 2 < Am 2 2 < 9.1 x 10" 5 eV 2 ; 0.23 < sin 2 6 U < 0.37 ; 

with the best fit for Am 2 2 = 8.1 x 10 -5 eV 2 and sin 2 612 = 0.29 (from solar and KamLAND 
data) . Note also that a previously acceptable reg ion with a somewhat higher mass splitting 



Am 2 2 ~ 1.4 x 10 4 eV 2 (the LMA II solution 30]) is excluded now at about 4a by the latest 
KamLAND data. 
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• finally, by using direct constraints from the CHOOZ reactor experiment as well as combined 
three-neutrino fitting of the atmospheric and solar oscillations, one can set the following 
upper limit on the #13 mixing angle: 

sin 2 9i 3 < 0.022 . 

The procedure we use in searching for a fit to neutrino sector parameters is as follows. First 
the low scale values of the quark and lepton masses and the CKM matrix angles and phase are 
chosen. (Generally we take a fixed value for mj and mt, while the other parameters are chosen 
randomly from a predefined range; however, mj and mt can also be chosen randomly). Next 
we pick a value for tan/3 and Msusy, and compute the quark- lepton sector quantities at GUT 
scale. Here we determine the relation between the lepton Yukawa couplings and quark Yukawa 
couplings, which amounts to determining the parameters a, b and phases a^fc, in Eq. ©. The 
phases combinations Z{ = a, + b{ are chosen as input (that is, they are picked randomly), while 
a, b, and the remaining two phases are obtained by the procedure of fitting the lepton eigenvalues 
described in Section III. Finally, we scan over the parameters which appear in the neutrino sector 
(if the neutrino mass matrix is either of type-I seesaw or of type-II seesaw, there is only one phase 
a; if both types appear, there will be two extra parameters, the relative magnitude and phase of 
the two contributions). 

The rest of this section is devoted to a detailed analysis of the predictions of the SO (10) minimal 
model for the neutrino sector, in type-I, type-II and mixed scenarios. (Due to its relative simplicity, 
we will start with the type-II case). However, let us first summarize our results. We find that in 
the type-II scenario, there is no good fit to the neutrino sector if the CKM phase is consistent with 
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experimental measurements (around 60 deg). This is in agreement with previous analysis 
however, our results are a bit more encouraging, in that that for 5c km = 7r /2 we find reasonably 
good fits, which improve significantly with not very large increases in the CKM phase. We can 
obtain marginal fits for Sqkm as low as 80 deg. More interesting are the results for the type-I 
case; here we can find good fits to the neutrino sector for values of 5c km as low as 50°, certainly 
consistent with experimental limits. As such fits have been not found before, one might consider 
this to be the main result of our paper. Also, we find that in the mixed case, there is possible 
to obtain a good neutrino sector fit in the case when the contributions coming from type-I and 
type-II are roughly equal in magnitude and of opposite phase. 
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A. Example of Type I Seesaw Fit 

We give here a representative example of a fit obtained in a type-I dominant case. This is 
obtained for m s {M z ) = 0.07 GeV, m c {M z ) = 0.85 GeV, m b (M z ) = 3 GeV, M t = 174 GeV, 
tan j3 = 40 and Msusy = 500 GeV. The values of the quark and lepton masses at GUT scale (in 
GeV) and the CKM angles are: 



m u 


= 0.0006745 


m c 


= 0.3308 


m t 


= 97.335 


m d 


= 0.0009726 


m s 


= 0.02167 


m h 


= 1.1475 


m e 


= 0.000344 




= 0.0726 


m T 


= 1.350 




= 0.2248 


•S23 


= 0.03278 


S13 


= 0.00216 






Sckm 


= 1.193 . 







(10) 



Here the masses are defined as m u = Y u smj3 vo, ra^ = l^cos/3 vo where Y u ,Y,i are the corre- 
sponding Yukawa couplings, and vq = 174 GeV is the SM Higgs vacuum expectation value 5 . The 
values of GUT scale phases (in radians) and a, b parameters are given by: 

a u =0.881 a c = 0.32678 a t = 3.0382 

b d = 3.63235 b s = 3.23784 b b = 0. (11) 
a = 0.08136 b = 5.9797 a = 3.244 . 

With these inputs, one can evaluate all mass matrices at GUT scale. In order to compute the 
neutrino mass matrix at M z scale, we use the running factors T22 = 1-06, = 1.03, where 



r22 



V M v33 )mJ\ J Mour ' V M, 33 J M z ' \ M.33 J . 



M z \ lvl v33/ Mgut \ lvl u2,3/ M z \ lvl u2>3/ Mqut 

with i,j = 1,2. The elements of the neutrino matrix above are evaluated in a basis where the 
lepton mass matrix is diagonal. 

One then obtains for the neutrino parameters at low scale: 

Aml 3 /Amj 2 ~ 24 , sin 2 9 12 ^ 0.27 , sin 2 26> 23 ^ 0.90 , sin 2 26>i 3 ~ 0.08 . 

Note here that only the atmospheric angle is close to the experimental limit, the solar angle and 
the mass spliting ratio being close to the preferred values. The elements of the diagonal neutrino 



5 One can write Eqs. (J^J, 1121 in terms of either the Yukawa couplings of the leptons and quarks, or their masses 
(that is, Yukawa couplings times running Higgs VEVs). In this paper we use the Yukawa couplings, but we multiply 
by the Higgs VEVs at the SUSY scale for simplicity of presentation. One can easily check then when going from 
one convention to the other, just the parameter a rescales, while b does not change. 
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mass matrix are 

m vi ~ {0.0021 exp(O.llt) , 0.0098 exp(-3.06i) , 0.048} 

in eV, with a normalization Am^ = 2.2 x 10 _3 eV 2 . The phases of the first two masses are 
the Majorana phases (in radians). Moreover, the Dirac phase appearing in the MNS matrix is 
4>d = — 0.23rad, and one evaluates the effective neutrino mass for the neutrinoless double beta 
decay process to be 

\'Y J XJl i m vi \ ~ 0.009 eV . 
B. Type II seesaw 

Much of the recent work on the neutrino sector in the minimal 5*0(10) has concentrated on the 
scenario when the type-II seesaw contribution to neutrino masses is dominant. The reason for the 
interest in this case is that, with: 

M v ~ M R ~ Mi- M d . 

b — r unification at the GUT scale, mj — m T , naturally leads to a small value of (M u )^ and hence 
large mixing in the 2-3 sector [6j]. However, while the general argument holds, it has been difficult 
(or impossible) to fit both large a ^ m and the hierarchy between the solar and atmospheric mass 
splittings at once. In this section we will try to show why this is so, and under which conditions 
this might be achievable. 

We will use the same conventions as in section III (that is, we work in a basis where M d 
is diagonal, and the parameters a, b and {M d )^ = nib are real and positive). However, in the 
construction of the neutrino mass matrices there will be an extra phase besides those which were 
relevant for the quark-lepton mass matrices. This phase a can be though as an overall phase of 
M[. One then has: 

M R = y{e i °M l - M d ) 
aM° = -{be ia + 2)Mie ia + ?>M d . (12) 

Following the analysis in Sec. Ill one can write: 

(M,) 2 2 ^ |6| m s e ib * 
(M,) 23 - a rm e W3 V 32 V33 

(M z ) 33 = am t e ia3 +bm b ~ m T e ia , (13) 
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with at close to ir and a + b s = e close to zero. Then, the neutrino mass matrix will be proportional 
to: 

. . / m s e i{e - a \b - e~ icT ) m 23 \ 

(M v ) [2i3] ~M x -M A d e-™ ~ ( 14 ) 

\ m23 m T e ta — m^e tcr J 

Note also that 77123 = m t e ia *V3 2 V33 is almost real positive, and due to the fact that amt — brrib 
and m(,V32 — m s , the 22 and 23 elements in the neutrino mass matrix are roughly of the same 
order of magnitude (in practice, one get 77123 somewhat larger than 77722). One then sees that if the 
phase a is chosen such that the two terms in the 33 mass matrix element cancel each other (that 
is a ~ —a) then there will be large mixing in the 2-3 sector, with: 

4. la \ m23 

tan(0 w )23 - • 

^33 

However, this is not the whole story. One needs also some hierarchy between the atmospheric 
and solar neutrino mass splittings: 

Am _k = r < 1 



AmL " ~ 20 



(based on the experimental measurements of neutrino parameters reviewed in the previous section). 
In terms of the eigenvalues 7772,7773 of the mass matrix (J14|) . one then has: 



m| ^ / |7772 2 m 3 3 -m| 3 1 \' < _L_ 

m 3 \l m 22l + l m 33l + 2|m 23 | 2 y ~ 20 ' V ' 

In order for this to hold, one needs a cancellation between the 7772277733 and 77723 terms in the 
numerator of the above fraction. This in turn imposes a constraint on the phases involved: 

(j) = Arg(777 r - m h e-^ a - a) ) ~ . (16) 

More detailed analysis shows that it is not possible (or very difficult) to get tan(^ I/ )23 larger 
than 1 while satisfying the relation (|15(1 between eigenvalues. However, this will create problems 
with the atmospheric mixing angle. The PMNS matrix is 

UpMNS = UiU u 

where Ui, U u are the matrices which diagonalize the lepton and neutrino mass matrices, respectively. 
Since the lepton mass matrix has a hierarchical form, the XJ\ matrix is close to unity, with (J7j)j 3 ~ 
(0, x, 1), where x = (77723/^33);*- The atmospheric mixing angle will then be: 
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FIG. 1: Graphical representation of the relationship between the b and r masses at the GUT scale and the 
magnitude and phase of the 33 element of the neutrino mass matrix. 

where the I and v lower indices make clear that we are discussing elements of the lepton and 
neutrino mass matrices. Note, however, that Eq. (|16|) implies that (77133)/ and (777,33);, have the 
same phase of m r ; then, since (77123)2 = ( m 23)i/j the net effect of the rotation coming from the 
lepton sector is to reduce the 2-3 mixing angle. Practically, since \x\ ~ 0.2, even if one has a value 
of tan(0^)23 close to one from the neutrino mass matrix, tan0 a ^ m will become of order 0.7 after 
rotation in the lepton sector is taken into account. 

This situation is represented graphically in Fig. As discussed above, eft = corresponds to 
the case most favorable for getting the right solar-atmospheric mass splitting ratio, while (f> = it 
corresponds to the case of maximal mixing angle. In practice this means that most reasonable fits 
are actually obtained when the angle </> is close to around 7r/2 (otherwise generally either the angle 
or the mass ratio are too small) 6 . Note however that <p = it (or any value greater than about 7r/2) 
would require that at GUT scale > m T . We may infer that in order to obtain a large mixing 
in the 2-3 sector one needs that at GUT scale m;, should be close to m T . 

This in turn can be insured by requiring large m\,{Mz) and/or tan/3. For example, in Fig. [2] we 
present the results obtained for values m\,{Mz) = 3.11 GeV, m t = 181GeV, Msusy = 500GeV and 
tan/3 = 50 (with these values, the ratio rv T = mb/m T (MGUT) — 0.96). Also, here we set the CKM 
phase 5 = 90deg, and let the other quark sector parameters vary between the limits discussed 
in section IIIA. The left panel shows the maximum atmospheric/solar mass splitting ratio R a / S 
as a function of the atmosperic mixing angle 9 a tm = #23- The three different lines correspond to 

6 Contributions from the phases in the lepton mass matrix can also improve the goodness of the fit (for example, if e 
is significantly different from zero, or a different from 7r). However, this generally requires that the parameters a, b 
have low values (as explained in section III). Hence we see that neutrino sector also prefers 8c km in the second 
quadrant and low m s - 
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FIG. 2: Left: maximum of the atmospheric/solar mass splitting ratio R a / S as a function of the atmospheric 
mixing angle 9 2 3, with cuts on the solar angle sin 2 26* 12 > 0.7 (dotted line), 0.65 (dashed) and no cut (solid). 
Right: maximum of R a / S as a function of the solar mixing angle #12, with cuts on the atmospheric angle 
sin 2 26> 23 > 0.9 (dotted line), 0.85 (dashed) and 0.8 (solid). 

different cuts on the solar mixing angle: sin 20i2 > 0.7 (dotted), sin 2 20 12 > 0.65 (dashed) and 
no cut (solid). One can observe here the correlation betewee large atmospheric mixing and small 
atmospheric-solar mass ratio. 

Conversely, the right panel shows the maximum of the ratio R a / S as a function of the solar 
mixing angle 9 so i = #i 2 . The three different lines correspond to cuts on the atmospheric mixing 
angle: sin 2 26> 23 > 0.9 (dotted), sin 2 26> 23 > 0.85 (dashed) and sin 2 26> 23 > 0.8 (solid line). We note 
here that the correlation between the solar angle and the mass ratio has the form of a step function 
(abrupt decrease in R a / S once sin 2 2#i 2 goes over a certain threshold), while there seems to be a 
close to linear correlation between the maximal solar and atmospheric angles. 

It is interesting to consider how these results change if the paramenters mb(Mz),m t , Msusy 
and/or tan/3 are modified. One finds out that the neutrino sector results have a strong dependence 
on the parameter tan (3. For example, if one keeps the parameters used in Fig. Infixed but increases 
tan /3, one finds that the fit for the atmospheric angle - atmospheric/solar mass ratio improves to a 
certain amount. That can be traced to the fact that the ratio r^/ T increases with tan/3. However, 
one also finds that the solar angle generally gets smaller. This happens because there is a correlation 
between the solar angle and the value of the s quark mass at GUT scale; namely #i 2 increases with 
r s j h = m s /mi J (MGUT)- On the other hand, the ratio r s / b decreases with increasing tan/3. 



17 




FIG. 3: Maximum of the ratio R a / S as a function of the atmospheric mixing angle #23 (left) and solar angle 
(right), for tan/3 = 40 (dotted line), 50 (solid) and 55 (dashed). Additional cuts are sin 2 26> 12 > 0.7 (left) 
and sin 2 26> 23 > 0.9 (right). 

Fig. exemplifies this behaviour. The three lines correspond to the maximum value for the 
mass splitting ratio R a i s , at values of tan (3 = 40 (dotted) , tan /3 = 50 (solid) and tan = 55 
(dashed line). A cut on the solar angle sin 2 20\2 > 0.7 is also imposed in the left panel, and a cut 
on the atmospheric angle sin 2 29\2 > 0.9 in the right panel. One can see that at larger values for 
tan/3 one might potentially get better fits for atmospheric angle and the atmospheric/solar mass 
ratio; however, the constraint on the solar angle becomes more restrictive. 

Smaller variations of the neutrino sector results will follow modification of the parameters 
mb(Mz),mt and Msusy- However, these variations follow the same pattern as above: that is, 
an improvement in the fit for the atmospheric angle due to the increase of the ratio r&/ r (which 
can be due to an increase in mt, or a decrease in Msusy) coincide with a worsening of the fit 
for the solar angle. As a consequence, the results presented in Figs. El ED can be improved only 
marginally. Scanning over a range of parameter space 2.9GeV < m^Mz) < 3.11GeV, 174GeV < 
m t < 181GeV, 500GeV < M SU sy < ITeV and 10 < tan/3 < 60 7 , we find the best fit to the 
neutrino sector to be sin 2 2^23 — 0.88, sin 2 2^12 — 0.74 and the atmospheric/solar mass splitting 
ratio R a / S — 24. We note that although these numbers provide a somewhat marginal fit to the 

7 In practice we find that the best results are obtained for large mt(Mz), large mt and large tan/3 (such that rt/ T 
is between 0.96 and 1). 



18 




FIG. 4: Maximum of the ratio R a / S as a function of the atmospheric mixing angle 6*23 (left) and solar angle 
(right), for S C km = 80° (dotted line), 90° (solid) and 100° (dashed). Additional cuts are sin 2 26 12 > 0.7 
(left) and sin 2 26» 23 > 0.9 (right). 

experimental results (the mixing angles are close to the exclusion limit, while the value for mass 
ratio is central) they are still allowed. 

However, the results discussed above were obtained for a value of the CKM phase 5 = 90deg 
which is too large compared with the measured value {5c km = 59 ± 13deg from PDG 

0)- As 

argued in section III (and indeed noted by previous analysis) there is a strong dependence on 
the goodness of the fit on the value of 5ckm, with larger values giving better fits. We show 
this dependence in Fig. The parameters are the same as in Fig |2] (m^Mz) = 3.11 GeV, 
rat = 181GeV, Msusy = 500 GeV and tan/3 = 50), but the three lines correspond to different 
values for Sckm' 80° (dotted line), 90° (solid) and 100° (dashed line). One can notice a rapid 
deterioration in the goodness of the fit with decreased 5ckm- Thus, for 5ckm = 80deg, the best 
fit to neutrino sector we find (after scanning over the SUSY parameter space) is sin 2 2023 — 0.88, 
sin 2 26»i2 ^ 0.7 and R a/s ~ 18. 

For purposes of illustration, we give a fit obtained for a type-II dominant case for 5ckm = 
80deg, m b {M z ) = 3.11 GeV, M t = 181 GeV, tan/3 = 55 and M SU sy = 1 TeV. The s,c quark 
masses at low scale are m s (Mz) = 0.074 GeV, m c (Mz) = 0.83 GeV. Then the values of the quark 
and lepton masses at GUT scale are (in GeV): 
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m u 


= 0.0008185 


m c 


= 0.3772 


ni t 


= 139.876 


m d 


= 0.0015588 


m s 


= 0.03554 


m b 


= 2.3547 


m e 


= 0.000525 




= 0.1107 


m T 


= 2.420 


S12 


= 0.225 




= 0.0297 


S13 


= 0.00384 






5cKM 


= 1.4 . 







The values of GUT scale phases (in radians) and a, b parameters are given by: 
a u = -0.4689 a c = -1.0869 a t = 3.0928 

b d = 2.6063 b s =2.2916 b b =0. (18) 
a = 0.09093 b = 4.423 a = 3.577 . 

The running factors for the neutrino mass matrix are r 2 2 = 1.09, T2z = 1.18. 
One then obtains for the neutrino parameters at low scale: 

Aml 3 /Am 2 12 ~ 18 , sin 2 20 12 ^ 0.7 , sin 2 26> 23 ~ 0.88 , sin 2 20 13 ~ 0.094 . 

The elements of the diagonal neutrino mass matrix (masses and Majorana phases) are 

m ui ~ {0.0016 exp(0.27i) , 0.011 exp(-2.86i) , 0.048} 

in eV. 

The Dirac phase appearing in the MNS matrix is 4>d = — 0.007rad, and one evaluates the 
effective neutrino mass for the neutrinoless double beta decay process to be 

\Y, U l m vi\ - °- 01 eV • 
C. Type— I seesaw 

The fact that in type-II seesaw one can obtain large mixing in the 23 sector is due to a lucky 
coincidence: the type-II neutrino mass matrix being written as a sum of two hierarchical matrices 
{Mi and M^), the most natural form for the neutrino mass matrix is also hierarchical. However, 
since 33 elements of both matrices Mi and Md are roughly of the same magnitude, by choosing the 
relative phase between the two to be close to ir, one can get a neutrino mass matrix of the form 
suited to explain large mixing in the 2-3 sector. 
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The question arises then if such a coincidence happens for the type-I seesaw neutrino mass 
matrix. To see that, let's write the Dirac neutrino mass matrix in the following form: 

rD be™ + 2 



Mi 



Mr + M d 



where Mr is the scaled right-handed neutrino mass matrix Mr = Mi — M d e ia and M d is a rescaled 
down- type quark diagonal mass matrix (the scaling factor in this later case £ = (b — e~ la ) / (be la + 2) 
is close to unity, since b is roughly of order 10). Then the type-I seesaw neutrino mass matrix 
would be: 

M vI = MfM^M? ~ M R + 2M d + M d M R ~ l M d . (19) 

Now, for most values of the phase a, Mr is hierarchical, therefore so is Mr 1 , therefore the 
type-I neutrino mass matrix is the sum of three hierarchical matrices (M d being diagonal). So it 
is not surprising that for most values of the phase a M u j is also hierarchical. What is remarkable 
is that there are some values of a for which the type-I seesaw mass matrix has a large mixing in 
the 2-3 sector, and moreover, this happens for the same values of a as in the case when the type-II 
mass matrix is non-hierarchical (that is, a close to ir). 

In order to see this let us consider the magnitude and the phase of the 33 elements (the largest 
ones) in the three terms on the right-hand side of Eq. (|19|). If a is not close to ir, the magnitude 
of MR33 is of order m&, with varying phase ((f) in Fig. Q in the first quadrant); the magnitude 
of (M d )s3 is also of order m&, and the phase ~ —a. For the last term, we make use of the fact 
that Mr being hierarchical, (Mr ^33 ~ 1/(Mr)ss ~ 1/W&; then M d M.R 1 M d ~ mb, with a phase 
close to —2a. We see then that for most values of a (M u j)^ is of order mj, while the off-diagonal 
elements are small. However, for a ~ ir, the cancellation in the 33 element of Mr is matched by 
a cancellation between the 33 elements of the 2M d and M d M.R 1 M d terms from Eq. (|19|) (since 
the relative phase between these is also a), thus leading to a non-hierarchical form for the type-I 
seesaw neutrino mass matrix. 

The fine-tuning between different contributions to the neutrino mass matrix is thus a little bit 
more involved in the type-I seesaw case compared to the type-II seesaw, but it can still lead to 
large mixing in the 2-3 sector. Moreover, since the correlations between the input parameters and 
the neutrino mass matrix elements are not so strong, most of the constraints discussed in the above 
section do not hold (for example, m& does not have to be necessarily very close to m T ). This may 
lead one to believe that it is possible to obtain a better fit for the neutrino sector in type-I models, 
and we found that in fact this is the case. 
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FIG. 5: Left: maximum of the atmospheric/solar mass splitting ratio R a / S as a function of the atmospheric 
mixing angle 6*23, with cuts on the solar angle sin 2 26> 12 > 0.8 (dotted line), 0.7 (dashed) and no cut (solid). 
Right: maximum of R a / S as a function of the solar mixing angle #12, with cuts on the atmospheric angle 
sin 2 26>2 3 > 0.9 (dotted line), 0.85 (dashed) and 0.8 (solid). 

For example, we show in FigJS^left) the maximum atmospheric/solar mass splitting ratio R a / S 
as a function of the atmosperic mixing angle Q a tm = #23 - with cuts on the solar mixing angle: 
sin 2 26*12 > 0.8 (dotted), sin 2 26\2 > 0.7 (dashed) and no cut (solid). In the left panel we show 
maximum R a u as a function of the solar angle for sin 2 2^23 > 0.9 (dotted), sin 2 26*23 > 0.85 
(dashed) and sin 2 2f9 23 > 0.8 (solid line). This fi gure is obtained for values m b {Mz) = 3.0 GeV, 
rat = 174GeV, Msusy = 500GeV and tan /3 = 40, while the CKM phase is allowed to vary between 
60 and 70 deg. We see that it is possible to obtain a large atmosperic/solar mass splitting ratio 
for values of the atmospheric and solar mixings consisted with experimental constraints. 

How do these results change if we modify the SUSY parameters tan/3 and Msusy, and/or 
the Mz scale masses? We find that one can get good fits for values of the parameter rj b / T = 
m b /m T (M GUT ) between 0.83 and 0.9. For values oim b {M z ) = 3.0 GeV, m t = 174 GeV, M SU sy = 
500 GeV, this means that tan/3 can vary between 10 and 55. However, if one increases m& or 
mt, generally rj(,/ T will increase (there is also a slight decrease with increasing Msusy, but less 
pronounced). Thus, for m b (Mz) = 3.11 GeV, m t = 174 GeV, one has to take tan/? roughly 
between 10 and 45 in order to get a good neutrino result. Fig. |H1 illustrates this behavior: the three 
lines correspond to values tan/3 = 20 (solid), 45 (dashed) and 55 (dotted), for m b {Mz) = 3.11 
GeV, m t = 174 GeV and Msusy = 500 GeV the corresponding values for are 0.86, 0.9 and 
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FIG. 6: Maximum of the ratio R a / S as a function of the atmospheric mixing angle #23 (left) and solar angle 
612 (right), for tan/3 = 20 (solid line), 45 (dashed) and 55 (dotted). Additional cuts are sin 2 26> 12 > 0.7 (left) 
and sin 2 26>23 > 0.9 (right). 

0.96. As noted above, the fit worsens dramatically for tan j3 > 45. Sokm is taken between 60 and 
70 degrees here also. 

Finally, the #13 mixing angle is found to lie in a range 0.06 ^ sin 2 2#i3 <■ 0.11, with a preferred 
value ~ 0.085 . We show in Fig. Qthe distribution of values for the #13 mixing angle and the Dirac 
phase in the neutrino mixing matrix 5n (in radians) obtained for a scan of parameter space with 
2.95 GeV < m b {M z ) < 3.05 GeV, 172 GeV < m t < 176 GeV, 500 GeV < M SU sy ^ 750 GeV, 
and 0.83 ^$ %/ r < 0.9 . The cuts on the other neutrino sector parameters are sin 2 2#23 > 0.88 (for 
all lines), sin 2 29 12 > 0.7 and R a / S > 18 (solid line), sin 2 26>i 2 > 0.75 and R a / S > 18 (dashed line), 
sin 2 2#i2 > 0.7 and R a / S > 20 (dotted line). Note that the preferred value for #13 is close to the 
experimental limit. 



D. Mixed type 

We end by considering the scenario when both types of seesaw are present, and give a non- 
negligible contribution. One can envision two different situations: first when one contribution 
dominates, and the other can be treated as a small perturbation, and second when the two contri- 
butions are of roughly the same order of magnitude. 

Let us first discuss the first case. It is obvious that the fits one obtains when considering 
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FIG. 7: Distribution of values for the #13 mixing angle (left) and the Dirac phase in the neutrino mixing 
matrix <5jy (right) obtained after a scan of the parameter space. The three lines correspond to different cuts 
on the other neutrino sector paramenters (see text for details). 

the 'pure' case (either type-I or type-II) can only improve. Indeed, by adding the other type 
contribution, we introduce another two free parameters: the relative magnitude and phase of the 
subdominant contribution relative to the dominant one. One can then obtain the results for the 
'pure' case by allowing the relative magnitude go to zero. 

Thus there appears an interesting question. Is it possible to improve the fit for the type- 
II scenario to a level which is compatible with experimental data? Unfortunately, the answer 
seems to be no. Our numerical simulations show that by adding type-I contribution as a small 
perturbation, one obtains a small improvement over the 'pure' type-II case, but not a significant 
one. 

This might seems counterintuitive at first glance. It would be reasonable to assume that if 
one takes the best result for a type-II fit and then adds another small quantity parametrized in 
terms of a free magnitude and phase, one should be able to obtain a somewhat better fit. The 
reason why this does not happen is that the type-I seesaw perturbation will modify in first order 
only the 33 neutrino matrix element 8 . However, as we saw from our analysis of type-II case, the 
goodness of the fit is determined mainly by the 22 and 23 matrix elements. That is, in type-II 

8 It is easy to see that in the sum of last two terms on the right-hand side of Eq. (1191 the 33 element is dominant if 
the type-II 23 mixing is close to maximal. 
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case, one generally has the freedom to easily modify the 33 element (by changing the phase a or 
the ratio m&/m T ) so in a sense for a best fit 77133 already has its optimal value, and adding a small 
perturbation to it does not buy any additional freedom. At most what one can achieve is to take 
some fit which is not optimal and make it better by adding the type-I perturbation. 

Let us now consider the case when the contributions from the two types of seesaw are of 
comparable magnitude. In this case, the neutrino mass matrix can be thought of as having the 
form displayed in Eq. (|19|). with the first term Mr (proportional to the type-II seesaw contribution) 
enhanced or diminished, depending on the relative phase between the two contributions. We find 
that if the type-II term is close to zero (that is, both contributions have about the same magnitude 
and the relative phase is close to 7r) then we can again get good fits for the neutrino sector. In 
particular, in this case we can obtain truly maximal mixing in the atmospheric sector (sin 2 2#23 — 1) 
with large enough mass splitting, and large enough solar angle. Also, obtaining these results does 
not require that we use a specific range of (low scale) parameters - like large tan j3, or large m&. 
Generally, as long as one can fit for the lepton-quark sector, one can obtain a fit for the neutrino 
sector, too. 

Here is a particular fit obtained in a mixed case. We took Sckm = 60deg, m^Mz) = 3.0 
GeV, M t = 181 GeV, tan/3 = 45 and Msusy = 500 GeV. The s,c quark masses at low scale are 
m s (Mz) = 0.071 GeV, m c (Mz) = 0.79 GeV. Then the values of the quark and lepton masses at 
GUT scale are (in GeV): 



m u = 0.0009440 m c = 0.3958 m t = 143.23 
m d = 0.001451 m s = 0.02991 m b = 1.6903 

m e = 0.00047 m M = 0.0992 m T = 1.910 (20) 

S12 = 0.220 S23 = 0.0320 s 13 = 0.00249 

Sckm = 1-05 . 

The values of GUT scale phases (in radians) and a, b parameters are given by: 
a u = -1.2527 a c = -0.39256 a t = 3.07385 

b d = 3.0896 b s = -3.1367 b b = 0. (21) 
a = 0.089727 b = 6.6365 a = 3.1505 . 
The running factors for the neutrino mass matrix are T22 = 1-04, riz = 1.09. We add the two 
contributions: 

M v ~ M uII j^Y + M vI re** , 
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with relative magnitude r = 0.9813 and phase <f> = 3.200rad. 
One then obtains for the neutrino parameters at low scale: 

Am| 3 /Am? 2 ~ 28 , sin 2 20 12 ^ 0.85 , sin 2 2023 - 0.98 , sin 2 26 13 ~ 0.05 . 

The elements of the diagonal neutrino mass matrix (masses in eV and Majorana phases) are 

m vi ~ {0.0028 exp(0.026i) , 0.0093 exp(-3.07i) , 0.0478} . 

The Dirac phase appearing in the MNS matrix is 4>d = — 0.3rad, and one evaluates the effective 
neutrino mass for the neutrinoless double beta decay process to be 

\J2 U ei m »i\ ~ °- 009 eV • 
V. CONCLUSIONS 

An 5O(10) model with only one 10 and one 126 Higgs multiplets coupling to fermions provides 
an appealing candidate for a unified theory at large scales. The number of free parameters in such 
a model is smaller than the number of free parameters in the Standard Model, thus giving the 
theory some predictive power. In particular, there appears the possibility that large mixing in the 
neutrino sector can be understood as a consequence of b — r unification at GUT scale. 

We revisit here the analysis of the minimal 5O(10) model and its implications for the neutrino 
sector. Our work differs from previous works in this area in several aspects. First, we consider the 
most general formulation of the model, with all the CP phases taken into consideration. Second, we 
use a new method for fitting the GUT scale parameters (Yukawa couplings) to the low scale masses 
and mixing angles in the quark and lepton sector. The running of the Yukawa couplings from low 
scale to unification scale is also taken into account, as well as dependence upon SUSY parameters 
like tan/3 and Msusy- We also analyze all possible cases for the neutrino mass generation, namely, 
when either the type-I or type-II seesaw mechanism dominates, or the case when both contributions 
are roughly of the same magnitude. 

Our results are as follows. For the type-II seesaw case, we find that the requirement for close 
to maximal mixing in the 2-3 sector, together with the large hierarchy between atmospheric and 
solar mass splittings, pushes the CKM phase to large values (reasonably good fits can be found 
for &ckm > 100°). This is in agreement with previous results. Interestingly, we also find the 
requirement that the solar mixing angle is large imposes significant constraints. Better fits can be 
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found at large tan/3, and large values of mb(Mz); however, the results obtained for values of Sckm 
in the first quadrant are at most marginally in agreement with experimental data. 

More interesting is the case when the type-I seesaw mechanism dominates. Here, contrary to 
previous analyses, we find that it is possible to obtain good fits for the neutrino sector for values 
of Sckm as low as 50°, which includes the range consistent with experimental results. Actually, 
unlike the type-II case, it seems that the goodness of the fit is not very dependent on Sckm] the 
relevant parameters seem to be the values of m c and m s at Mz scale, with preference for larger 
values for m c and lower values for m s . The solar angle also imposes significant constraints on 
the parameter space, but in this case it is possible to obtain good fits for a larger range of tan/3, 
m b {M z ) and M t . 

It is also interesting to note that, for values of the phase a appearing in the neutrino sector close 
to 7T, the type-I and type-II neutrino mass matrices are roughly proportional. This means that if 
the type-I mass matrix is non-hierarchical (thus leading to large mixing in the 2-3 sector), so will 
the type-II mass matrix be. Hence, if one considers the case when both contributions from type-I 
and type-II must be taken into account, one sees that the dominant contribution will determine 
the type of the fit. However, a special case is when the two contributions are roughly of the same 
magnitude, and their relative phase is close to it. In this case they cancel each other, and the 
resulting neutrino mass matrix is generally non-hierarchical. One can obtain very good fits for the 
neutrino sector in this case, too. 
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